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ABSTRACT 

We are interested in black holes in Loop Quantum Gravity (LQG). We study the simple model of 
static black holes: the horizon is made of a given number of identical elementary surfaces and these 
small surfaces all behaves as a spin-s system accordingly to LQG. The chosen spin-s defines the area 
unit or area resolution, which the observer uses to probe the space(time) geometry. For s = 1/2, 
we are actually dealing with the qubit model, where the horizon is made of a certain number of 
qubits. In this context, we compute the black hole entropy and show that the factor in front of 
the logarithmic correction to the entropy formula is independent of the unit s. We also compute 
the entanglement between parts of the horizon. We show that these correlations between parts of 
the horizon are directly responsible for the asymptotic logarithmic corrections. This leads us to 
speculate on a relation between the evaporation process and the entanglement between a pair of 
qubits and the rest of the horizon. Finally, we introduce a concept of renormalisation of areas in 
LQG. 
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I. INTRODUCTION 

Loop Quantum Gravity (LQG) is a canonical quantization of General Relativity, which relies on a 3+1 decomposition 
of space-time (for reviews, check lj). It describes the states of 3d geometry and their evolution in time (through the 
implementation of a Hamiltonian constraint). The states of the canonical hypersurface are the well-known spin 
networks, which represent polymeric excitations of the gravitational field. The main and first achievement of the LQG 
framework is the implementation of quantum area operators and the derivation of their discrete spectrum. A surface 
is now made of elementary patches of finite quantized area. Based on this discrete structure of the 3d space, one can 
study in details the entropy associated to a surface and the entanglement between the patches making the surface. The 
main goal of the present work is to study analytically the relation between entropy and entanglement of the horizon 
state and discuss their physical/geometrical interpretations. 

Let us discuss more precisely the concepts that were mentioned above. Entanglement can be loosely defined as an 
exhibition of stronger-than-classical correlations between the subsystems. For a long time, it stood up among the 
apparent "paradoxes" of quantum mechanics for its nonlocal connotations Q . Recently it became one of the main 
resources of quantum information theory 0,03. We discuss its relevant properties in Sec. IIB. In our work we are not 
concerned with its role in quantum information processing. Instead we are interested in its relations to the entropy of 
black holes and the role of entanglement as a pointer to their physical properties, similarly to the recent discussion of 
quantum phase transitions 

A spin network is a graph, a network of points with links representing the relations between points. Each link or 
edge is labeled by a half-integer, which defines the area of a surface intersecting the link. More generally, a surface 
is defined through its intersections with the edges of the spin network described the underlying quantum 3d space: 
the surface can be thought as made from elementary patches, each corresponding to a single intersection with a link 
and whose area is given by the label of that link. These labels actually stand for SU(2) representations i.e spins, and 
are conventionally noted j. Due to regularization ambiguities, there is not a definitive consensus on the precise area 
associated to a given spin j. The original and conventional prescription is of an elementary area a(j) = \J j(j + 1), 
but other reasonable choices are a(j) = j or a(j) = j + 1/2 (although one would like to keep a(j = 0) = in the end 
for both physical intuitive reasons and also mathematical consistency [jj). Now, points, or vertices, of a spin network 
represents chunks of volume. Mathematically, they are attached a SU(2) intertwiner - a SU(2) invariant tensor between 
the representations attached to all the edges linked to the considered vertex. 

A generic surface on a spin network background is thus described as a set of patches, each punctured by a unique link 
of the spin network. A SU(2) representation j is attached to each patch. We denote its Hilbert space V 3 . Intuitively, 
a (quantum) vector \jm) in the Hilbert space corresponds to the geometrical normal vector to the surface defined by 
the elementary patch. Now the spin network defines how the patches, and therefore the whole surface, is embedded 
in the surrounding 3d space and describes how the surface folds. For a closed surface, the region of the spin network 
which is inside the surface exactly defines an intertwiner, invariant under SU(2), between the patches of the surface. 
This intertwiner contains the rotation-invariant information on the way the patches organize themselves to form the 
surface e.g the angle between two patches. 

An important remark is that any spin j representation V 3 can be decomposed as a symmetrized tensor product of 2j 
spin-i representations V 1 ^ 2 . Therefore, one can interpret that a fundamental patch or elementary surface is a spin-^ 
representation. All higher spin patches can be constructed from such elementary patches. For example, considering 
two spin-i patches, they can form a spin representation or a spin 1 representation: in one case, the two patches are 
folded on one another and cancel each other, while in the later case they add coherently to form a bigger patch of spin 
1. Considering an arbitrary surface, one can then look at it at the elementary (fundamental) level decomposing it into 
spin-i patches, or one can look at it at a coarse-grained level decomposing the same surface into bigger patches of spin 
s > 1/2. From this point of view, the size of the patches used to study a surface is like the choice of a ruler of fixed size 
used by the observer to analyze the properties of the object. When studying the entropy and entanglement on a given 
surface, we will therefore consider it made of a certain number of patches 1/2 at the elementary level (and not allow 
the size of the patches to vary when performing entropy computations). Then one can study the coarse-graining or 
renormalisation of these quantities when one observes the surface at a bigger scale, using bigger patches to characterize 
the surface. 

In the present work, we would like to study some basic properties of Black Holes in the framework of Loop Quantum 
Gravity. Black holes are actually the main object studied in quantum gravity. Indeed, on one hand, they are especially 
simple objects from the point of view of theoretical general relativity, and on the other hand we expect a consistent 
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theory of quantum gravity to shed light on the origin of the entropy of a black hole and on the so-called "information- 
loss paradox" [|| Q • Here we propose to analyze the basic features attached to the horizon of a black hole through a 
simple model inspired from Loop Quantum Gravity. 

We deal with the simplest case of a Schwarzschild (non-rotating) black hole at the kinematical level. Actually 
studying black holes at the dynamical level in quantum gravity should be very interesting, and should allow to 
understand the evaporation process at the level of space-time. We will later comment on the dynamics as induced by 
LQG, but this will not be the main point of our study; instead we have decided to focus on the black hole kinematics. 
Considering a static black hole, our main assumption is that the only information accessible to the observer outside 
the black hole can be seen on its (event) horizon. Considering the horizon as a closed surface, LQG describes it as 
made of patches of quantized area and describes the interior of the black hole in terms of (a superposition of) spin 
networks whose boundary puncturing the horizon defines the patches of the horizon surface. From the point of the 
external observer, it does not matter what is inside the black hole, but only the information which could be read 
off the horizon is relevant. More mathematically, the horizon surface is defined as a set of elementary patches to 
each of which is attached a given SU(2) representation. Let us label the patches with an index i = 1, ..,n and note 
the corresponding spins jt. The space geometry within the horizon, i.e. the black hole, is fully described by a spin 
network, which can possibly have support on a complicated graph within the horizon. From the point of view of the 
horizon, the details of the spin network inside the black hole do not matter: what actually matters is which intcrtwincr 
between the SU(2) representations on the horizon does the spin network induce. Indeed the observer outside the black 
hole is only interested by the state of the patches on the boundary, which lives in the tensor product V^ 1 ® .. (g> V^ n . 
The only constraint on the state is that it should be globally gauge invariant, i.e. SU(2) invariant, such that the 
horizon state is exactly described by an intertwiner between the representations . From this point of view, horizon 
states are described by a spin network having support on a graph with a single vertex inside the horizon, i.e. a fully 
coarse-grained spin network. Actually, we would like to stress that our analysis applies to generic closed surface states 
as soon as we decide to ignore the details of the geometry of the enclosed spatial region. 

The basic model for a quantum black hole that we propose to study is not new. We assume that the model describe 
the black hole horizon at the fundamental level. In other words, we assume that the observer describing the black hole 
have the maximal possible resolution and probes space(time) with a ruler of minimal area i.e. corresponding to a spin 
1/2 surface. The resulting model for the horizon is to consider it made of a number n of spin 1/2 patches: a horizon 
state will live in the tensor product (y 1 / 2 )®". More exactly, we want a SU(2) invariant state, i.e. an intertwiner 
between these n representations. At this point, let us point out that a spin 1/2 system is usually called a qubit (in the 
language of quantum information Q), so that one can say that the horizon is made of n qubits and that a horizon 
state is a singlet state for these n qubits. Finally, the area of the horizon surface will by definition be n x ai/2 where 
ai/2 is the area associated to an elementary spin 1/2 patch. 

In this model, we will present the calculation of the entropy of the black hole and recover the entropy law, in an 
asymptotic limit, with a first term proportional to the horizon area and a logarithmic correction with a factor —3/2. 
For these purpose, we develop a method to compute the number of intertwiners as a random walk modified with a 
mirror at the origin. Then the point that we would like to particularly stress is that the n qubits on the horizon are 
correlated, and more precisely entangled: the horizon is not made of n uncorrelated qubits and the —3/2 factor is 
actually a reflection of the fact that these qubits are entangled. 

Then one can renormalize, or equivalently coarse-grain, this fundamental model, now assuming that the observer 
probes the space(time) with a spin s > 1 elementary surface. A coarse-grained model of the black hole is then to 
consider horizon states to be given by intertwiners between n representations of a fixed spin s. One can repeat the 
entropy calculation and check how the entropy law gets renormalized. The first term proportional to the area should 
simply get scaled by the ratio of the area a s corresponding to a spin s surface to the ai/2- Then we will prove using 
the random walk analogy that the logarithmic correction factor is universally —3/2 and does not get renormalized. 

Our approach is different from the standard loop quantum gravity approach, which studies the classically induced 
boundary theory on the black hole horizon (or generically any isolated horizon). It is shown that one gets a U(l) 
Chern-Simons theory which one can quantize and count the quantum degrees of freedom. In particular, some 
ways of computing the black hole entropy in such a context also involve parallels with random walk calculations 0. 
Nevertheless, one would ultimately like to start with the full quantum gravity theory, identify horizons within the 
quantum states of geometry and study the quantum induced boundary theory on the identified quantum horizon. 
However, searching for horizons in a spin network states requires either introducing a well-defined notion of observer 
at the quantum level and/or being able to extract the (semi-classical) metric from the quantum state of geometry. We 
do not tackle these issues in the present work. We decide to model the black hole as a region of the quantum space 
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with a well-defined boundary (the horizon) such that the only information about the geometry of the internal region 
accessible to the external observer is information which can be measured on the boundary. Our results can therefore 
be applied to any closed surface in loop quantum gravity. Issues about actually identifying this closed surface as a 
horizon at the semi-classical level is left for future investigation. 

The paper is organized as follows. The next section will detail the fundamental spin-i black hole model. We will 
compute the entropy of the black hole and entanglement between qubits on the horizon. This will lead us to relate the 
entanglement of a pair of qubits (or more generally a small part of the horizon) with the rest of the horizon with its 
probability to evaporate from the black hole. We remark that the evaporation process can not work on a single qubit 
so that minimally only pairs of spins 1/2 can evaporate. Also, introducing a time scale canonically associated to a 
black hole, like the (minimal) time of flight of a particle on a circular orbit around the black hole, we derive the rate 
of evaporation from the probability computed previously. Section III deals with the generalized spin s models. Using 
the random walk analogy, we derive the entropy law and show the universality of the 3/2 factor for the logarithmic 
correction. Section IV tackles the issue of the renormalisation of surface areas in the framework of LQG. Consider a 
surface defined microscopically on a spin network. The spin network describes how the surface is folded and embedded 
in space. A natural issue is what will be the coarse-grained area that an observer will assign macroscopically to that 
surface. We use the calculations of intertwiners to deduce the most probable area seen macroscopically assuming no 
knowledge of the underlying spin network. We then propose to refine this by defining a (background independent) 
state of a surface in (loop) quantum gravity through the probability amplitude assigned to intertwiners between the 
elementary patches making the surface. This information can be equally thought as data on the surrounding spin 
network or as data on how the surface geometry gets coarse-grained (how its geometry looks at different scales). 



II. THE SPIN-1/2 BLACK HOLE MODEL 

We describe the horizon of the black hole as a two-sphere punctured by the underlying spin network by 2n, «eN, 
edges labeled by spins i. The number of punctures defines the area of the horizon: A = 2na 1 / 2 , where a 1 / 2 is the 
area corresponding to a spin-i elementary patch. The geometry of the interior of the black hole is described by the 
potentially complex graph inside the two-sphere. However, this information is not available to an observer living 
outside the horizon. This external observer has only access to the information on the horizon, i.e., to the horizon 
state. It is given by the intertwiner between the 2n punctures on the horizon, which corresponds to completely coarse- 
graining the internal spin network into a single vertex. The state lives in (y 1 / 2 )® 2 ™ anc l i s SU(2)-gauge invariant. The 
tensor product of 2n spins | simply reads: 

2n n 

0^ = 7 i 8ffiii . (II .i) 

3=0 

where a n ,j is the irreducible representation of the permutation group of 2n elements corresponding to the partition 
[n + j,n — j]. More details of the representation theory of the permutation group can be found in Appendix D. The 
dimensions = dim(cr nj ) count the degeneracy of each spin j in the tensor product (y 1 / 2 )® 2 ™. The intertwiner 
space H is defined as the SU(2)-invariant subspace of the tensor product i.e. the spin j — subspace. It is obviously 

isomorphic to a n ,o- We note its dimension N = Let us remind that we are considering an even number of spins 
since there does not exist any intertwiner between an odd number of 1/2-spins. 

It is known that the degeneracy factors are given in terms of the binomial coefficients: 



In particular the dimension of the intertwiner space is 



N=--—Cl l . (II.3) 

n+1 2,1 y 1 



Below we give a derivation of these formulas in terms of random walks, which will be straightforwardly generalizable 
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to higher spin models when we replace the fundamental 1/2-spin by an arbitrary s-spin. Using the Stirling formula 1 
to compute the asymptotics of the Catalan number C^ n l '(n + 1), we obtain the asymptotic behavior of the entropy of 
the quantum black hole: 

S = In N ~ (2n) In 2 - - ln(rc) - - In n. (II.4) 

n~ >oo 2 2 

The first term is linear in n and thus in the area, and is the usual term that one expects to appear in any black hole 
entropy calculation. The precise factor in front of the area A actually depends on the unit used to define areas. The 
next term is the logarithmic correction in ln(„4). It is independent of the unit used to measure areas and thus the 
—3/2 factor should be universal (invariant under coarse-graining). We will not discuss the other terms, which can be 
considered as negligible in the large area case. In the next section, we will show that we get the same factor in any 
s-spin black hole model. 

Under the assumption that the external observer has no knowledge of the internal geometry of the black hole, the 
density matrix ascribed to the black hole horizon is the totally mixed state on the intertwiner space: 

1 N 

p=nY,\ i ^> ( IL5 ) 

i=l 

where form an arbitrary orthonormal basis of the intertwiner space. It is obvious that the von Neumann entropy 
of this state reproduces the (micro-canonical) entropy: S(p) = — trplnp = InN. Below we study the entanglement 
properties of this state and show how the logarithmic correction to the entropy formula reflects the correlations between 
the 2n patches making the black hole horizon. 

Let us end this introduction by discussing the generality of this 1/2-spin model. Thinking about rotating black 
holes, some might be tempted to only require invariance of the state under U(l) rotations around a particular axis 
-say J z - and thus count all states with vanishing angular momentum m = instead of the stronger constraint of 
vanishing spin j = 0. Then the dimension of the horizon state would be C£ n and we would have only a —1/2 factor 
in front of the logarithmic correction. However the SU(2) invariance which we require is the gauge invariance present 
in loop quantum gravity and has nothing to do the physical SO(3) isometry group of the two-sphere. A more serious 
criticism would be that in loop quantum gravity one must allow spin network labeled by arbitrary spins and thus we 
should more generally allow punctures of arbitrary spin and not restricting them to 1/2-spins in order to have the true 
picture of a quantum black hole. Our counter-argument is that any spin k can be decomposed into 1/2-spins and thus 
allowing higher spin punctures would lead to an overcounting of the intertwiners and horizon states 2 . More precisely, 
for example, an intertwiner between one 1-spin and (2n — 2) 1/2-spins automatically defines a unique intertwiner 
between (2n) 1/2-spins (through the unique intertwiner V 1 / 2 ® V 1 !" 1 — > U 1 ). From this perspective, one can consider 
the 1/2-spin black hole model as the universal black hole description in loop quantum gravity. We will discuss this issue 
in more details in the later Sec. lIVI on area renormalisation in loop quantum gravity. Let us nevertheless point out that 
the interested reader can find details on the structure of the intertwiner space when allowing arbitrary spins labeling 
the punctures in [l2T |. Another argument in favor of the 1/2-spin black hole model is given in p"H as a symmetry 
argument in the semi-classical limit. This is related to the fact that the intertwiner space provides a representation of 
the permutation group, which can be interpreted as the discrete diffeomorphisms of the black hole horizon as made 
as 2n elementary patches. We will tackle this issue later in Sec. III. 2. 



We remind that Stirling formula gives the asymptotical behavior of n\: 




2 One can argue that replacing two 1/2-spins by one 1-spin actually modifies the area. This depends on the precise microscopic area 
spectrum. If one chooses the obvious linear spectrum cij oc j, then one avoids this issue. Nevertheless, the standard prescription in loop 
quantum gravity is a,j oc \J j(j + 1) which gives 

01 = Af 
2oi /2 V 3* 

If this case, the coarse-graining actually modifies the area. One should then keep in mind that a rigorous entropy count should involve 
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II. 1. Entropy from random walk 



In this section, we show how one can compute the degeneracies of the tensor product decomposition and thus the 
entropy as a random walk calculation. Indeed, describing the decomposition of (y 1 / 2 )® 2 ™ iteratively, the analogy 
with the random walk is clear. If we are at the /c-spin representation V k , tensoring it with V 1 ^ 2 , we obtain the 
representations of spin k — 1/2 and k + 1/2. Only if we are at the trivial representation k = 0, then V° ® V 1 / 2 is 
simply V 1 ! 2 . This way, it appears that the degeneracy of the fc-spin in the tensor product (y 1 / 2 )® 2 ™ can b e computed 
as the number of returns to the spot k after 2n iterations of the random walk with a mirror (or a wall) at the origin 
0. More precisely: 

Proposition 1 The degeneracy coefficients in the decomposition of the tensor product (i/ 1 / 2 )® 2 ™ { n to £/j e irre- 
ducible spin-j representations equal the number of returns RWM(j) to the spot j after 2n iterations of the random 
walk of step 1/2 with a mirror in j = and starting at the origin. Moreover the random walk with a mirror is given 
in terms of the standard walk by the simple relation: 

RWM„(j) - RW„(i) - RW„(j + 1), (II.6) 

where RW(j) is the number of returns to the spot j after 2n iterations of the random walk. 

To prove this result, it is more convenient to deal with integer steps instead of switching between half-integers and 
integers after each step. Then we actually consider the tensor product as (("i/ 1 / 2 )® 2 )®™; a ^ eacn step, we are tensoring 
with (y 1 / 2 )® 2 = V° V 1 i.e. moving either one step up or one step down or staying at the same spot. Then it is 
straightforward to check that RWM„(> j) = J2k>j RWM„(fc) satisfies the same iteration law as the number of returns 
RW n (j) of the standard walk without mirror at the origin. 

Then as RW„(j) = , we recover the previous formula dj — — C^^ +1 . In particular, this shows the 

simple identity: 

J24 n) = RW„(0) = C 2 "„. (II.7) 

j 

Having these exact formulas, it is straightforward to derive the asymptotics of RW„(j = 0) and RWM„(j = 0) when 
n grows large, and the result is 2 2 ™/ v / n and 2 2n /n^/n, respectively. 

On the other hand, it is possible to directly extract these asymptotics without computing the exact expression. 
Indeed, working on the unit circle U(l) and using the orthogonality of the modes e m0 , one has: 

RW„(j) = — f +7T d9e- t(2 ^ e (e ie + e- l9 ) 2n = 2 2n x — / + " dO e - <(2j)e (cos 0) 2n . (11.8) 

For j — 0, we evaluate this integral in the asymptotic limit n — > oo using the saddle point approximation. First, we 
write 

RW„(0) = 2 2n x - f +2 d9 (cos 6) 2n = 2 2n x - [ + * dO e 2n log(cos 9) . (11.9) 

Then the exponent (f)(9) = log(cos 9) is always negative and has a unique fixed point (which is actually a maximum) 
at 9 = with the second derivative (fj^(9 = 0) = — 1. Therefore we can approximate the integral by: 

i r+i 2 9 2n 1 r + °° 2 

RW„(0) cx 2 2n - / d9e~ ne ~_x-/ dxe~ x , (11.10) 



counting all the states of a given classical area A up to an uncertainty 8A. The dependence of the chosen scaling of 5A with A would 
then certainly change the details of the entropy law. 
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and derive the right asymptotic behavior. Following the same line of thought, starting from 

RWM n (j = 0) = RW„(0) - RW„(1) = 2 2n ±- f + dd (1 - e- ae )(co S 9) 2n , (11.11) 



and using that this expression is real, we get: 



RWM„(0) = 2 2 ™- / 2 d6 (1 -cos 26») (cos 9) 2n = 2 2n - f * dd (1 - cos26»)e 2nlog(cose) . 

Now we have an integral of the type J <Wf(ff)e~ n ^ e \ The exponent <f>(9) has only one fixed point at 9 = 0. Then 
since f{9) vanishes at 9 — 0, the leading order in the asymptotics is given by the first non-vanishing derivative of /: 



2 



+ 00 



RWM„(0) oc 2 Zn - / d99 A e- rM — - / dxx l e~ x . (11.12) 

n—>oo TT J n^' 2 7T J _ QO 

This shows the 2 2n jr?l 2 asymptotic behavior of the dimension of the intertwiner space and thus of the horizon state 
space. The same analysis will be generalized in the next section to any s-spin horizon model. The factor —3/2 of 
the logarithmic correction to the area-entropy law will appear to be universal and, as shown above, has a simple 
explanation in terms of random walks. 

Although the asymptotics in 1/y/n and X/riy/n have a natural interpretation from the random walk perspective, 
a simpler way to obtain integral formulas for the degeneracies is to use the orthogonality of the SU(2) characters. 
Indeed, we have: 



d?° = 



dgxM{Xi(9)) =- / sin 2 9d9x j (9){xi(0)) , (H.13) 



where dg is the normalized Haar measure on SU(2). The character Xj{9) IS the trace of the group element g in the 
j-spin representation and is a function of only the rotation angle 9: 

, x sin(2? + 1)6* 
xm= V g , X|(0)=2cos0. 



II. 2. Entanglement 

The black hole model we analyze allows to deal with entanglement in a finite-dimensional setting. In addition, we 
consider only a bipartite entanglement, i.e., entanglement between two distinguishable parts of this system. This is 
the best understood part of the entanglement theory 0, Il5j| . In this section and in the appendices we present only the 
most essential elements of the analysis. The missing proofs and the generalizations to other families of states can be 
found in |l6j . 

Pure state entanglement is easily identified. Pure states that cannot be written as a direct product of two states 
are entangled. For example, the singlet state of two qubits, 

l*-> = ^(l + ->-!-+», (H.14) 

is an entangled state. The definition is more involved for mixed states. It can be expressed as a part of the threefold 
hierarchy. At the lowest level there are direct product states p = pi ® p2- Then, there are classically correlated (or 
separable) states, 

p = ^w l p\®p\, ^u>i = l, Vw 4 > 0. (11.15) 

i i 

which are mixtures (convex combinations) of the direct product states. Finally, entangled states are those can not 
be represented as separable. They have their own hierarchy which does not concern us here. Apart from two-qubit 
systems, no single universal criterion that allows to tell whether a mixed state is entangled or not is known. 
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For pure states there is a natural way to quantify entanglement. The degree of entanglement of a pure state \^S) is 
the von Neumann entropy of either of its reduced density matrices, 

E{\y)(y\) = S(p AtB ) = -tTpA,Blog P A,B, p A = tr B \*)(% (11.16) 

For example, reduced density matrices of |\& _ ) are the maximally mixed spin-| states pa.b — l2xa/2, so its degree of 
entanglement is log 2 = 1 bit. In general, a pure state whose degree of entanglement is logd, where d is the dimension 
of the Hilbert space of one of the subsystems, is called maximally entangled. 

There are various measures of entanglement of mixed states that reflect different aspects of their preparation. We 
adapt here the entanglement of formation, which is defined as follows. For all possible decompositions of the state 
Pab = P into mixtures of pure states, 

p = J2 w i\^i)(^il X>* = 1 ' Vwi>0, (11.17) 

i i 

the weighted average of degrees of entanglement of the constituents is calculated, and the minimum is taken over all 
decompositions 

E F (p) = mmY2wiS(pi). (11.18) 

For pure states this expression reduces to the degree of entanglement. Entanglement of formation is zero for unentan- 
gled states. For generic states, an analytic expression for Ep(p) exists only for two-qubit systems. 

Nevertheless, we were able to calculate entanglement of p of Eq. (|II.5(1 for all bipartite splittings (Proposition 2 
below). Moreover, it is possible to show that all standard entanglement measures coincide on p [l6[, so there is no 
arbitrariness the choice of the entanglement measure. Henceforth we drop the subscript and simply write E(p). 

Let us split the 2n qubits into a group of 2k and a group of 2n — 2k qubits and adapt the angular momentum basis. 
Hence the Hilbert space is (y 1 / 2 )® 2k ^ ^yi/2)®(2n-2fe) ) for which a shorthand H = H a ®'Hb will be used. The basis 
states in either of the subspaces are labeled as \j, to, aj) and \ j, to, bj), respectively. Here < j < k and —j<m< j 

A B 

have their usual meaning and aj enumerate different irreducible subspaces V 3 . Their multiplicities d- ' are given in 
Appendix B and aj — 1, . . . , d^ and bj — 1, . . .d?. The constraint J 2 = ensures the states |Xj) are the singlets on 
Va <£> V B subspaces for j — 0, .., k, 



1 



3 



\j, a 3 ,b 3 ) ee - 1 == ^ (-iy- m \j, -m, a,) \j, to, b 3 ). (11.19) 



-j 



A more transparent notation is based on representing each subspace with a fixed j as V J ® D J , where V J is a 
2j + 1 dimensional space that carries a spin-j irreducible representation (irrep) of SU(2) and D° is the degeneracy 
subspace. In the case of qubits the Shur's duality identifies these spaces as irreps of the permutation group S^n 
that correspond to the partition [n + j,n — j] of 2n objects [UGH. Hence each basis state can be represented as a 
tensor product state on three subspaces, 

\j,aj,bj) = \j) ab ® \aj) D i A ® Ibj)^ (11.20) 

In this notation Eq. I|II.5(I becomes 

P=^EE li>%') & i)0>i> & j'l = \3)tiUB®\aj){aj\ D A (8) \bj)(bj\ D B (11.21) 

j a,j,bj j aj,bj 

Reduced density matrices of these states, 

Pj.a 3 =tr B \j,aj,bj){j,a,j,bj\, (11.22) 

where here and in the following the tracing out is over Hb, are independent of bj. There are exactly d^ of the matrices 
Pj.a-j i so the reduced density matrix pa — tr bP is given by 

^ = ^EE d iW ( IL23 ) 

j a 3 

With this notation, it is possible to prove the following result: 
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Proposition 2 The entanglement of the black hole state p is given by 

E (p) = 4 E E ) = i E d f # J°g(jy + 1)- (n.24) 



The proof can be found in Appendix A. Following [T^ |. this result can actually be generalized to more general zero-spin 
states: 

Proposition 3 Consider an arbitrary convex combination of the zero spin states, 

P= T,jY, aj ,b j w ^fi j \3)^\AB®\as){a j \ J ^®\b J ){b s \ D a, 



£^^=1 (H-25) 
all measures of entanglement for the state p are equal to 

E(P)= E w ^ lo g(2j + l). (H.26) 

As an example we consider the entanglement between the block of 2 qubits and the rest. Dividing the system into 
2 and 2n — 2 qubits leads to the following degeneracies: 

i = 1, d» =( 2n ~ 2 Y-, df= ( 2n ) * n -V -, (11.27) 
011 \n-l)n 1 \n J 2(2n- l)(n+ 1) v 7 

that asymptotically satisfy <fjf /N ~ 1/4, df /AT ~ 3/4. Hence 

£(p|2) = ^|log3~ Jlog3. (11.28) 
Entanglement of 4 and 2n — 4 qubits is hardly more challenging, with asymptotic values of the degeneracies being 



As a result 



4 4) = d£d$/N~ 1/8, 4 4) =dfd£/N~9/16, s ( 2 4) = rf^df /JV - 5/16. (11.29) 

£7(p|4)~ Alog3+^log5. (11.30) 
16 16 

Fig. 1 illustrates the convergence of the entanglement to its asymptotic value. 

One finds that if the In qubits of the horizon are separated into two sets of n qubits each, then the entanglement of 
the state p (asymptotically) equals to E(p\n) ~ \ logn, while the entropy S(p) the reduced density matrices p n of its 
halves is S(p n ) ~ nlog2. More technically, this translates to the quantum mutual information [17l between the black 
hole and its halves being three times the entanglement between the two halves: 



I p {n : n) = S{p n ) + S(p n ) - S(p) ~ 3E(p\n). (11.31) 



Let us recall that for any pure entangled state like the singlet state of Eq. I|II.14|) . we have: 

I*{A : B) = 2J5(*). (11.32) 

In a mixed entangled state the classical correlations are present on top of the entanglement, so it is natural to expect 
that I p {A : B) > 2E(p\A : B). From the quantum information-theoretical point of view, it can be argued that the 
quantum mutual information represents the total amount of classical and quantum correlations 01 ■ Hence our result 
shows that the logarithmic correction to the black hole entropy represents the total amount of correlations for the 
symmetric splitting of the horizon state. 

More generally, while the exact entanglement between the two parts of the spin network depends on k, an explicit 
calculation shows that Eq. i|11.31fl can be extended to 

I p {A : B) w 3E(p\A : B) (11.33) 

for any nearly arbitrary partitions into two parts, A and B. For example, for k = 10 the coefficient is approximately 
2.956. 
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FIG. 1: Entanglement between 2k and 2n — 2k qubits for 2k — 2 (left) and 2k — 4 (right). The entanglement of the two qubit 
block for n = 4 is 0.706 ebits, and the entanglement of the four qubit block for n = 6 is 0.9429 ebits. 



II. 3. About evaporation 

Let us now focus on the possible physical meaning of the entanglement calculations computed above. Let us consider 
a pair of qubit. It is in a mixture of having an intertwiner j — or j — 1 with the rest of the horizon patches. In 
the case that it has an intertwiner j = 0, standard loop quantum gravity tells us that the pair of qubits is actually 
not linked -detached- from the rest of the horizon qubits i.e. it is not on the horizon anymore. From the information 
point of view, a j = intertwiner means that the pair of qubits is uncorrelated with the rest of the horizon, so 
intuitively it can not be part of the black hole horizon. The entanglement, or more precisely the entangled fraction, 
is a precise quantification of how much the pair of qubits is in the intertwiner state j = 1 compared to being in in 
the intertwiner state j — 0. It is therefore natural to assume a relation between the unentangled fraction and the 
evaporation probability of the pair of qubits (more precisely that there exists a monotonous function relating the two). 

Let us consider the unentangled fraction of a pair of qubits s 2 ^ . It depends of course on the total number of patches 
n, i.e. on the area of the black hole horizon. 

The probability of evaporation of the pair of qubit is given by the Born rule: 

Probevap = Tr(fiP j=0 ), (11.34) 

where p is the black hole density matrix and Pj=o the projector on the subspace Hj=o of intcrtwiners having a link 
of spin j = between the chosen pair of qubits and the rest of the horizon. Then p is actually the identity on the 
intertwiner space H normalized by the dimension N = dimH, so that the probability is exactly computed as the 
unentangled fraction: 

rl^Tn-W ,A=2qubits >B=rest , , 

Prob - dlm7 ^° - d " d " - s ( 2 ) (II 35) 

Prob cvap - dimH - N -s . (11.35) 

Then to go from probabilities to evaporation rate we need a time scale. In our simple model, we do not have any 
dynamics. The only way to introduce a time scale is to go to the semi-classical limit: for large black holes, the natural 
time scale is proportional to its mass raMa vC4 oc y/n. For example, we can consider the black hole as a clock 
with its period being the time of flight of the shortest stable circular trajectory around the black hole. Finally, we can 
conclude that: 

d-4 An 4 2) 4 2) 4 2) 

— — oc — cx — — oc oc — (11.36) 

dt dt t M y/n K 1 



where the proportionality coefficients actually depend on the precise choice of time scale and the exact relation between 
n, the area A and the mass M. Since s converges to a fixed non-zero value 1/4 when n grows to infinity, we get 
a (approximately) fixed probability of evaporation (of the pair of qubits) for large black holes. To make the analysis 
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(2) 

more exact, we should consider the exact dependence of the uncntangled fraction s Q on the black hole size n. Actually 
Eq. (|II.27f> gives us: 

J 2 ) 



which leads to a correction to the usual evaporation formula: 



dA dn 1 3 

(X — — oc — - 



dt dt 4^ 8nv^' 

This might allow us to derive precise corrections to the evaporation formula of black holes. Nevertheless, we should 
then also consider the corrections to the semi-classical mass formula M tx ^/n. Moreover all relies on the use of a 
semi-classical time scale, which would need to be derived from the exact dynamics of the underlying quantum gravity 
theory At the end of the day, we recover in our simple model the standard Hawking evaporation formula: 

dM 1 
~df K ~M2- 

Our analysis is actually very similar to the one performed by Bekenstein and Mukhanov |Tflj . 

We would like to point out that we have made the implicit assumption that the evaporation process is completely 
dominated by the evaporation of a pair of qubits. Indeed, to be thorough, we should consider as well the possible 
evaporation of bigger blocks. Nevertheless, we have computed the unentangled fraction of a block of 2k qubits (we 
only consider blocks of even size because of the SU(2) invariance requirement): 



,(2fc) 



1 



k(n — k) 



And it is clear that the unentangled fraction drops rapidly with the block size k. Actually, only one out of 2 2k basis 
states of the 2&-qubit block does not contain smaller unentangled pieces, hence the "proper" unentangled fraction 
(probability of a block of size k to get detached without any of its smaller piece getting detached itself) is 2~ 2k s^ k \ 
which now drops exponentially with the block size k. 

We would like to insist on the fact that our analysis is not dynamical: we rely on the Born rule and treat the black 
hole evaporation as a traditional decay (like for radioactive emission). To go further, one should consider the geometry 
state of the interior and exterior of the quantum black hole, and not only the horizon state, and analyze their quantum 
dynamics. Our main point here is that it does make sense to relate the notion of entanglement between parts of the 
horizon to probabilities of evaporation from the horizon. 

A last point is that the pair of qubits which evaporates will likely not be the Hawking radiation: it forms a part 
of spacetime which is in the exterior of the black hole. The Hawking radiation will be produced by the difference 
of energy between the system of a black hole with horizon size 2n and the system of a black hole with horizon size 
2{n — 1) plus the detached pair of qubits. 

Finally, we have related the entanglement calculations between parts of the horizon to the Hawking evaporation 
of black holes. As entanglement is a purely quantum notion, this fits with the fact that the evaporation is a purely 
quantum effect. Moreover, it points toward a link between (quantum) information concepts such as correlation and 
entanglement and the physics of (quantum) geometry. 



II. 4. Black holes at the interplay between quantum gravity and quantum information 

As we have illustrated above, quantum information concepts, such as entanglement, are useful for the study of 
black holes in quantu m g ravity. This relation goes beyond conceptual points and borrowing the technical tools of the 
entanglement theory |la |. In particular, the use of SU(2) invariant subspaces in quantum black hole has a precise 
analogy in the establishing of quantum communication without a shared reference frame .20] . 

Communication without a shared reference frame (SRF) can be formulated as follows. Assuming that there are 
no sources of noise, different orientations of the reference frames of the communicating parties induce a unitary map 
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between their states and operators that describe their experimental procedures. This unitary representation depends 
on the physical nature of the qubits, usually massive spin-1/2 particles. When the parties have no knowledge about 
the respective orientation of their frames, the state of the transmitted qubits, the initial state p looks averaged over 
SU(2) as far as the receiver is concerned: 

T(p)= f dUUpU\ (11.38) 

JSV(2) 

where dU is the Haar measure of SU(2). The task is to transfer information in the most efficient way despite this 
decoherence. 

Indeed, the transmittable information is to be encoded in a SU(2)-invariant way. For this purpose, to start with, we 
can use the j = subspaces of the space of the n transmitted qubits. For example, in the decomposition of (V 1 / 2 )® 4 , 
there are two distinct j = subspaces, which can be be used to encode the states |+) and |— ) of a single qubit. In 
the notation of Eq. (|II.19(I it is given, e.g., by 

|+>-»|i = 0,1,1), |-)-»|j = 1,1,1). (11.39) 
As we have seen earlier, the Schur's duality |13| gives in general: 

n n/2 

W„ = (g)C 2 = V j ®C7 n j. (11.40) 

A SU(2) transformation on H. n will act on the spaces V 3 while leaving the degeneracy spaces <7 nj invariant. It is then 
straightforward to check that the map T acting on density matrices on TL n will randomizes the V 3 sector (mapping 
to the maximally mixed state) while preserving the information stored in <J n ,j- More precisely, it is possible to show 
the following: 

Proposition 4 Consider the space H. n and the map T acting on the space of density matrices overTL n - Let us consider 
the space H.2n = <S> Hn an d its S\] (2) -invariant subspace H^n = °~2n,j=o- Then Im(T) is given exactly by the 
reduced density matrices defined as the partial traces over B of density matrices on Ti^n ■' 

Tib (Den(W°j) = Im(T„). 

In particular, for p € Im(T„), there exists a p density matrix on Ti^n such that p = Trs(p). Then by definition, the 
entropy of p is S(p) = E(p\A : B) the entanglement of the state p. 

The most efficient transmission method is a block encoding, when L logical qubits are encoded into n physical ones. 
To this end the qubits are encoded into the degeneracy labels of the highest multiplicity spin-j max representation of 
SU(2). As shown in Sec. 4, the highest multiplicity spin is asymptotically j max = \/n/2. The number of logical 
qubits is the logarithm log 2 of the dimension of the degeneracy space corresponding to j max , hence the efficiency L/n 
asymptotically approaches 1 as 1 — n^ 1 log 2 n. 

The same setting can be applied to cryptographic purposes |2lj . Imagine that Soraya tries to communicate to 
Marc and that they share a common reference frame to measure and send qubits. On the other hand, we consider 
an eavesdropper Tamara who does not have any information about that common reference frame. In this context, for 
all state p that Soraya sends, Tamara will see T(p). One can define the quantum capacity of the secured quantum 
channel following [21j . We are looking for the biggest possible Hilbert space 7i C 7i n such that for all density matrix p 
with support in 7i, T(p) is the maximally mixed state po on 7i and Tamara will not be able to extract any information 
about the message which Soraya wants to send to Marc. The quantum capacity is then, as usual, the logarithm of the 
dimension of the maximal such Hilbert space TL. 

This is easily done using the tools introduced previously. If po is the maximally mixed state on TL then log(dim7Y) = 
S(po). Moreover, there exists a state po such that po = TybPq, and then the entropy is given by the entanglement 
S(po) = E(po\A : B)). Finally, proposition Ogives us: 

n n 

^) = EE<W 2 -? + 1 )' with EE<k- = L 

j—0 aj ,bj j—0 cij ,bj 
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Obviously, the maximal entanglement will be when w^ n ' = 1 and = for all j < n: Soraya must use the highest 
possible spin j = n to encode her messages. Then the quantum capacity is given by the entanglement, so that 
QCap = log(2n + 1), as obtained in plj . 

Finally, [2l| notices that the classical capacity of this secured quantum channel is three times its quantum capacity. 
This factor 3 is exactly the one between the classical correlations and quantum correlations which we obtained earlier 
in equation l|II.33|) . 

To conclude this section, we have shown that this particular protocol of quantum communication without shared 
reference frame allows to translate directly physical properties of the quantum black hole, like its entropy or entan- 
glement, to quantum information notions such as classical and quantum channel capacities. We hope that such a 
dictionary could be further developed. 

III. THE SPIN-1 BLACK HOLE MODEL AND GENERALIZATION 

We now describe higher spin models, where we assume that the horizon is made of n elementary patches of a fixed 
spin s. We start by analyzing the 1-spin case and then we generalize the result about the entropy of the general 
case. These models can be interpreted as coarse-grained models of the fundamental 1/2-spin model when the observer 
counting the degrees of freedom does not have access to a resolution finer than the unit area a s set by the elementary 
patch of spin s. 

In the spin-1 case, we are interested into the decomposition of multiple tensor products of V 1 : 

n n 

&V 1 = ® 1 df>Vl. (IH.1) 

3=0 

We show below how to compute the degeneracies l dj in terms of random walk. Then we explain how the random 
walk calculation extends to the general case of the decomposition of multiple tensor products of V s : 

n n 
3=0 

This allows us to prove that the black hole entropy S^ s \n) = log( s d n ^) always have a first term proportional to 
the horizon area and then a logarithmic correction with the same —3/2 factor. Similarly, Eq. 1|II.33|) that compares 
quantum mutual information I p (A : B) for an arbitrary bipartite splitting of the horizon spin network with the 
entanglement between its halves holds also in this case. 

This claim was already made in [Toj |. However, there, the authors start with the number of conformal blocks for the 
q-deformed SU(2). We believe that our approach is simpler and more transparent, and that the random walk analogy 
allows a straightforward generalization to a large class of gauge groups (other than SU(2)). For example, the same 
results with the universal —3/2 factor directly applies to the supersymmetric extension osp(l|2), which is the relevant 
group for the loop quantization of N = 1 4d supergravity. 

III.l. Random walk analogy, entropy renormalisation and universal log correction 

Following the same logic as for the 1/2-spin model, we describe the decomposition of the tensor product <g> n V s 
iteratively. Since for an arbitrary spin j representation we have 

V s ®V 30 = v j , 
i=lio-s| 

the iterative process can be thought as a random walk. From the spot jo, one can go in one step with equal probability 
to any spot from jo — s to jo + s, as long as jo is larger than the fixed spin s. When jo is smaller than s, then we run 
into the wall at the origin j = 0. We formalize this in the following statement: 
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Proposition 5 The degeneracy of the irreducible spin representation Vj in the decomposition of the tensor product 
<g) n V s is the number of returns to the spot j after n iterations of a random walk allowing jumps of length up to s with 
a wall/mirror at j = 0. This truncated random walk is such that the number of returns can be simply computed from 
a standard random walk allowing jumps of length up to s with no obstacle at j = 0: 

s df = RWM« (j) = HW« (j) HW« (j + 1). (III.2) 
Then the number of returns of the standard random walk can be computed as an integral: 

RW 1 S) W = ^ J +7! d0e- ije (e- lse + .. + 1 + .. + e +ls6 ) n , (III.3) 

for s is an integer and 

RW« (j) = d9 e-^ e (e-^ e + e -*(2*-2)* + .. + e -w + e « + .. + e -H(2.)fl) , ( ln .4) 

wften s is a half-integer. 

As in the 1/2-spin case, it is actually possible to derive these integral formulas directly using the characters of SU(2): 



s df = J dg X M(Xs(9)) n = I £d6 sin9 sin(2j + 1)6 ( g^|±^ ", (HI.5) 
For s = 1, we get explicit formulas of RW*(0) and RWM n (0) as sums: 



where [j/J is the integer part of j/. Distinguishing the cases with n is odd or even, one can further simply the expression 
of the truncated random walk 3 . We can directly evaluate the asymptotics of these numbers of return using Stirling 
formula. Indeed for large n, letting x = 2l/n € [0, 1], the sums can be approximated by the integrals: 

1 f , < „Af^ ™™ ,i /„n 1 f , 3a; — 2 n 4>(x) 

e i 



HWj,(0) ~ -L / i e^, RWMi(0) ~ J- / - 

where the exponent is given by 

4>{x) = a; In 2 — xh\x — (1 — x) ln(l — a;). 



It is straightforward to check that <f> admits a unique maximum at x = 2/3 for which 0(2/3) = In 3 and <j>"(2/3) = 
—9/2 < 0. We then directly extract the asymptotical behavior of the integrals using the saddle point approximation. 
Since l/xy/l — x doesn't vanish at x = 2/3 we get: 

on 

RW*(0) ex — , (III.8) 



3 For example when n = 2m + 1, we have: 
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while since (3x — 2)/a; 2 v'l — x vanishes at x = 2/3 but has a non-vanishing second derivative 4 , we get: 

on 

RWt<(0)cx — =. (III.9) 
n^Jn 

A faster and easier calculation is using the saddle point approximation directly on the integral expression given in 
the previous proposition [S] Indeed: 

1 r+w i /" + if ( i\n + ? 

RWi(0) = — / d6(l + 2cOs6) n = — / d9 e ™ ln (2co S e+l) + V J_ / ^ e nln(2cos0-l)_ (ill.io) 

2vr,/_ 7r 2ttJ-^ 2vr ./_* 

Both exponents (2 cos 6* + 1) and (2cos# — 1) have a unique maximum at 9 = 0. Thus the first term will contribute 
3 n /y/ri to the asymptotic behavior while the second term will only give a term of the order of 1/ ' \/n. In the end, we 
simply recover RW„(0) ~ 3 n /y / n. Now, 

RWM I 1 1 (0) = — / + d9(l -e- l6 )(l + 2 cos 6>)™ = — / + d9 (1 - cos 0)(1 + 2 cos 6) n . (III.ll) 

2?T J_ 7r 27T 

As (1 — cos 6*) ~ 9 2 /2, the saddle point approximation directly yields RWM„(0) ~ 3 n /n^/n. 
This proof can be straightforwardly adapted to the generic spin-s case, so that we get: 

Proposition 6 We have the following asymptotical behavior of the degeneracy of the trivial representation V° in the 
decomposition of the tensor product ® n V a : 



RWi*'(0) oc (2s t_ 1)n , s 4 n) =RWMW(0) oc (2s+ ^ )? ' . (111.12) 

n— >oo nWn 



Since the degeneracy of the trivial representation is simply the dimension of the intertwiner space i.e. of the Hubert 
space of horizon states, we derive the entropy law: 

S { n s) ~nln(2A+l)-|lnn + .., (111.13) 

where we recognize the first term proportional to the horizon area and the now universal —3/2 factor in front of the 
logarithmic correction. 

Here, we showed that the 3/2 correction factor comes directly from the analogy with the random walk and does 
not depend on the particular model (choice of spin for the black hole elementary surfaces) that we chose. This proves 
that the 3/2 factor is actually invariant under coarse-graining i.e under the choice of area unit that we use to probe 
the horizon of the black hole: this is therefore a universal factor and rigid prediction from these particular models of 
quantum black holes. 

Moreover, we can refine our evaluation of the degeneracies 8 dj pushing further the approximations of the integrals 
giving RWM: 

, An) 3 ^3~ (l + 2a)" 4s(l + s)n-3-6j(l+j) (1 + 2s)"(l + j) 

j j<n 8 V 2tt [s(l + s)p/2 n 5 / 2 1 +J) n~oo (s(l + s)n) 3 / 2 ' 

This allows us to show the following statement. 



4 Let us consider the integral f^ dx f(x)e n( t'( x h Assuming that <j> has a unique maximum at xq £]0, 1[, we write <f>(x) = a — b(x — xq) 2 
with b > 0. Then we expand / around xq and doing the change of variable to y = (x — xo)y/n, which gives the asymptotic behavior of 
the integral as: 

„na f + oo / „. 1 7 .2 \ „ 

E — / dy (f(x ) + f'(x )^L + U"(x ) 1 L + ..)e~ b y 2 . 

It is clear that when /(xq) vanishes, the first derivative }'{xq) doesn't play into role in the approximation since / dyye~ b y = 0, so that 
the first relevant term is the second derivative f"(xo). 
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Proposition 7 We have the following asymptotical behavior of the maximal dimensionality s dj — maxj s dj of 
the representation spaces V J in the decomposition of the tensor product ® n V s : 

j max oc Vn, s df ] oc (2S + 1) " . (111.14) 
The proportionality coefficients depend on s. 

The calculations leading to these results can be found in Appendix B. 



III. 2. Elementary scale and irreducibility 

The structure of spacetimes as given by LQG is fundamentally discrete. Hence the diffcomorphism symmetry that 
is defined on a smooth manifold is replaced by a discrete symmetry on a spin network ^l|. The Hilbert space of LQG 
is spanned by the spin networks, and it carries a representation of the permutation group. Any representation of the 
permutation group of n objects S n can be either faithful or irreducible. Moreover, for large n the squared dimension 
of (C fc )®" is smaller than n!, so the faithful representation of the permutation group on the spin network any size is 
impossible. There is no a priori reason |llj| to require that all the edges of a spin network are labeled by the same 
spin s, since it is not necessary for building representations of the permutation group. However, we can see that the 
requirement of the irreducibility of the representation is equivalent to describing spin networks only in terms of the 
fundamental representation s = |. Namely [l3L |22|. the Shur's duality that allows to decompose (C N ) n as a direct 
sum of tensor products of the irreps of SU(iV) and Sk for the direct product of U(iV) holds only when the dimension 
of the elementary space is N. For SU(2) it takes the form 

2n n 

(g)C 2 ^0^®(7 B|i , (III. 15) 

3=0 

where V J is the irreducible spin-j representation of SU(2), and the irreps a n j correspond to the partition [n+j,n — j] 
of 2n objects There is no similar structure for the direct product of different vector spaces on the left hand side, 
and the isomorphism in the case of the direct product of involves the irreducible representations of SU(iV). 



IV. ON AREA RENORMALISATION IN LQG 

Let us look at a generic surface, made of 2n elementary patches of spin- 1/2. Unlike the previous case we allow 
for open surfaces. Then the Hilbert space of surface states is simply the tensor product (g) 2 ™"!/ 1 / 2 . The precise state 
describes how the geometry of the surface, the way it is folded. It is the surrounding spin network, in which the surface 
is embedded, which provides the information on the surface state and the way that the surface is precisely embedded 
in the 3d space. 

Now the surface can be folded any way at the microscopic scale. Then depending on how the surface is folded, an 
observer at a larger will see a smooth surface which is only a approximation of the "real" surface and will measure an 
area smaller than the full microscopic area .Amicro = 2na\/2- 

Let us start by the case when we completely coarse-grain the surface state. That is we decompose ® 2n V x l 2 = 
®™dj , and the total spin j will describe the macroscopic area of the completely coarse-grained surface. For 
example, for 2n = 2, V x l 2 ® V 1 ^ 2 = Vq © V 1 and we have two coarse-grained area states: j = which corresponds to 
the surface folded on itself and j = 1 which corresponds to the open surface with the two elementary patches side by 
side. 

For large n, one can describe the completely coarse-grained state for the trivial state of geometry i.e. when we do 
have no information at all the state of the surface. Then the most probable macroscopic area is given by the spin j 
with the largest degeneracy dj . Given the explicit expression of these degeneracy factors as dj = C^,^ — C^^ +X , 
we easily write that 

4^ (rc-j)(2j + 3) 

d (n) ( n+i+2 )(2 i + l)' 
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and identify the maximal entropy for: 

(IV.l) 

One can reproduce the same result by using the Stirling formula approximation 5 to dj for large n. This square-root 
law is natural from the point of view of the random walk analogy and it is actually true in any s-spin model, when 
looking at the maximal degeneracy in the decomposition of the tensor product ® 2n V s (see Appendix B for details). 

However, as soon as we have a non-trivial geometry state and have more information of the tensor product state 
corresponding to the surface at the microscopic level, we do not expect the macroscopic area to follow this rule 
anymore. Nevertheless, we would like to insist on the fact that the maximal area given by the microscopic area 2nai/ 2 
is totally improbable at the macroscopic level: its degeneracy is simply 1 while typical degeneracies grow as 2 2n /n^/n 
and the maximal degeneracy as 2 2n /n. Therefore we believe that the study of area reorganization in LQG will not be 
as straightforward as assuming a simple linear rescaling of the microscopic area. 

More generally, we do not need to go directly at the completely coarse-grained level. Physically, we are interested 
by the whole coarse-graining process: for example going from the fundamental scale set the 1/2-spin to a larger semi- 
classical scale given a fixed large s-spin. A simple calculation one can do considering 2n qubits is to take a subset of 
m qubits and compute the reduced density matrix of the subset from the state (possibly mixed) of the full surface. 
Then one will know the most likely area of this subsurface formed by the considered 2k qubits. Pushing this set-up 
further, we can partition the full surface into p bigger patches formed by n/p 1/2-spins and compute the reduced 
density matrices of each bigger patch and thus obtain a coarse-grained description of the surface. A first step would 
be to do this calculation in the framework of the totally mixed state corresponding to the black hole state. 

Indeed, a reduced density matrix for a 2k qubit patch is 

k d t 

pm = jf y Ei, d fps*r ( IV - 3 ) 

3=0 a.j=l 

At this stage of the analysis it is not clear which prescription to calculate a coarse-grained value of j should be used. 
There are three reasonable options that all lead to jcoarsc ~ Vk, while the numerical factors are of the order of unity. 

• The expectation value of the patch's spin is 

(J) =trJ P(2k) , (IV A) 

where J = VlP or any other suitable function of J. Since all the density matrices are diagonal in the angular 
momentum basis, this expectation value can be calculated as 

<^> = E J 0>?' fc) > Pf' k) = df k) df n - 2k) /N, (IV.5) 

3 

and, e. g., J — ■J j(j + 1). We calculate the expectation value of J by replacing the sum in Eq. I|IV.5|) by the 
integral. Using the asymptotic expansion of Eq. l|IV.2(l the normalization of the probability distribution p 1 "™'^ 




Letting x = j/n S [0, 1], we have 

d\ n > « d {n \x) = -—f(x)e- n ^ x \ (IV. 2) 

3 \rnm 



with 

1 2x+ - 

f(x)= - — f^-, <p(x) = (l-x)ln(l-x) + (l + x)ln(l+x). 

Vl-x 2 1 + x + ± 

The maximum is given by the equation rufi' (x) = (Inf)'(x). For a; ma x w 0, we get Xmax ~ 1/V 2n. Moreover if we prefer to look at the 
maximum of (2j + l)df ] , then one simply multiply the function / by 2x + l/n and then the maximum is located at x max ~ 1/i/n. 
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is fixed by 



N', 



where fj, = k/(n — k) « fc/n, all the constants where absorbed in N and X max 
J ~ j the expectation value becomes 



(IV.6) 

\j\p2k — > was used. Taking 



N' 

Since for large spins the elementary area is approximately proportional to the spin, 



A 11 



as illustrated on Fig. 2 below 
A 
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500 1000 1500 2000 

FIG. 2: The macroscopic area of the 2n = 2000 qubit horizon as a function of the number of patches, p = 2, 4, 5, 
square-root area spectrum is assumed, a,- = \fj\j + 1), and (J) is taken to represent the patch. 



,2000. The 



Considering (J) as an expectation values of a spin measurement, the most likely spin value j that occurs in such 
a measurement is jmax that maximizes the probability distribution Pj . When 1 -C k the spin of the highest 
degeneracy representation is given by Eq. I|C.7|) . and for k <C n it reduces to 

= Vk. (IV .9) 



J 



max 



• The most likely state \j,j z ) (again summing over all the degeneracy labels) corresponds to j that maximizes 

(n.k) 

& = (iv-io) 

The value of j can be determined similarly to j max , and for 1 -C k <C n it is 

J=y/k/2. (IV.ll) 

For all the three choices of the coarse-grained j the area of the patch is 

.4(2*0 = «i/2«Vfc, (IV.12) 
where a ~ 1 depends on the coarse-graining procedure. As a result, the coarse-grained area of the surface is 

Anacro = pA( 2 k) = 0.i/ 2 a^/np. (IV. 13) 

In the end, it would be interesting to generalize these area renormalisation calculations to more generic situations in 
loop quantum gravity, and thereby underline the distinction between the microscopic geometry and the semi-classical 
measurable and observable geometry. 
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V. 



CONCLUSIONS AND OUTLOOK 



We have studied simple models for black holes in the framework of loop quantum gravity. We consider the black 
hole horizon as a quantum surface made of a certain number n of identical elementary surfaces. These elementary 
quantum surfaces are mathematically modeled as spin-s systems, while the whole horizon must be a singlet state 
(or intertwiner) in the tensor product of these n spins s. Actually our framework applies to generic quantum closed 
surfaces as soon as we coarse-grain the interior geometry. 

We have computed the entropy of black holes in this context and recovered the usual entropy with a logarithmic 
correction. The factor in front of that correction is actually —3/2 whatever the spin s which we choose. We also 
analyzed the entanglement between parts of the horizon. More technically, we computed the entanglement for splitting 
of the horizon into two parts of arbitrary sizes. We have shown that this entanglement is responsible for the logarithmic 
correction to the entropy law. More precisely, the logarithmic correction is the total correlation between the two parts 
of the horizon and is equal to three times the entanglement. Moreover, we also point out a potential relation with the 
evaporation process, which allows us to recover at first order the Hawking evaporation formula. From our perspective, 
quantum black holes appear at the interplay between quantum gravity and quantum information, and we show how all 
the objects appearing in our analysis of the black hole can be translated to the framework of quantum communication 
without shared reference frame. 

Finally, we have used all the developed technology to illustrate a concept of renormalisation of areas in loop quantum 
gravity: given a surface at the microscopic level, it might be folded in ways not observable to a macroscopic observer, 
so that the macroscopically observed area is different to the microscopic area postulated by loop quantum gravity. 
And we formulate a square-root law for the renormalisation flow of the area. We see a few possible development of 
the present work: 

• One can try to introduce the dynamics into these simple black hole models, using the quantization of the ADM 
energy in loop quantum gravity. This would allow to study the evolution of the entanglement, potentially look 
at the entanglement decay and the evaporation process. 

• One can try to take into account the rotation of black holes in our models. However, this means to introduce more 
objects, relative to which we would define the rotation of the black hole. This would need to model the exterior 
of the black holes and think about spin network geometries whose semi-classical regimes are approximately the 
Schwarzschild metric. 

• It would be interesting to push further our analysis of area renormalisation, and study this concept in more 
complicated situation that the black hole state which turned out to be rather simple. We could also look at the 
renormalisation of other geometrical quantities such as the volume. 

• Finally, we can push further the use of entanglement and correlations as tools to probe the geometry of spin 
networks in LQG. One can naturally speculate on a link between the notion of distance and the actual correlations 
between subsystems of a g iven spin networks, similarly to what happens in spin lattices. Work on this topic will 
be reported elsewhere [24j . 



We would like to thanks Laurent Freidel, Lee Smolin and Florian Girelli for their interest into the present work. We 
are also grateful to Steve Bartlett for a very very long discussion on reference-frame-based cryptography protocols in 
quantum communication during some warm days in Sydney. 



Here we prove the Proposition [21 giving the exact formula for the entanglement of the black hole state. In the basis 
we adopted the states \j, ctj,bj) are given be the bi-orthogonal (Schmidt) decompositions Q, so their reduced density 
matrices are diagonal. Moreover, the only nonzero part of pj >aj is a sequence of 2j + 1 terms l/(2j + 1). Using the 
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decomposition of A into and D* subspaces, 

pj, aj = ^p^[ l vi ® \aj)( a j\Di (A.l) 

For the future convenience we introduce pj = lyj/(2j + 1). Two matrices pj, a and pi^ ai have orthogonal supports if 
one or both of their indices are different. The reduced density matrix pa = tr bP is 

P a = nJ2J1 d fPi l a i><°i|j* • (A.2) 

Since the von Neumann entropy of the reduced density matrices pj satisfies 

S(pj,«,) = 5(p J ) = log(2i + l), (A.3) 
the weighted average of the entanglement in this decomposition is 

A B 

m c j c j m 

Mp) ^EEE ^) = ^E d X lo s(2j + l). (A.4) 

3=0^=11^ = 1 3=0 

All pure states \*f/ a ) that appear in alternative decompositions of p ought to be some linear combinations of the states 
\j,aj,bj), 

|*o> = c «,iijh\j)AB ® \aj) D f A ® \*>j) D 3 B - (A.5) 
The diagonal form of p forces the coefficients c a j aj bj to satisfy the normalization condition 

^w a c a ^ a]b] c* alaih = ■^5 j l6a l a j 8b j b r (A.6) 
a 

The reduced density matrices pa{o) = tr b| x I'q)(5 , q | are 

Pa{o) = Yl ''"-J 1 ' b l'< ® \ a i)( a 'j\D A - (A-7) 

3 i'. •',.«' 

Introducing 

^ajajttj = Ca ^ a i b i C a,ja' j bj l ^3 Q-Q = ^ (A.8) 

we rewrite the reduced density matrix of | v E' Q )( v I / a| as 

PA («) = J] ^ ® A J («) . ( A-9) 

3 

where the fictitious density matrix A was introduced on the degeneracy subspace D 3 A , 

A ^ a ) = Yj A «.^i< K>KU ■ ( A - 10 ) 

J a.a. 

From the orthogonality of the matrices pj ja , and Eq. (|A.6|I it is easy to see that 

E^^j","; = N d f s <*3°-' s - ( A - n ) 

a 

The weighted average of the entanglement of the decomposition {$ Q } is (5({$ Q })) = J2 a w aS{pA{o))- From 
Eq. (jA.lOfl and the concavity property of entropy 0, ^ follows that 
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<S({$ Q })) > PfoO + ^A^a))] > Y,*Oc<*j(a)S( Pj ) = ^^dfdfSfa). (A.12) 



Hence Se{p) is indeed the entanglement of formation. The equality of all measures of entanglement for the black hole 
states is proven in |16| . 

APPENDIX B: EVALUATION OF THE DEGENERACIES s ^"' FOR ARBITRARY SPIN 



We consider the integer spin. The half-integer spin is treated in exactly the same way. From Eq. (|III.3(I it follows 
that 

RWM;(;0 = i £ d6f(j,6)e n ^ F ^ e \ (B.l) 
where f(j, 9) = cos(jO) — cos[(j + 1)9} and 



2 cos | 

F(s, 9) = — (B \ ' > - 1. (B.2) 



Hence 



For a fixed j and n — > oo the saddle point method gives the following estimate of the degeneracy: 

sAn) „ 3 [Y (1 + 2sT 4.(1 + s)n - 3 - 6j(l + j) 

s 4 n) - ( ;^ 2s) T^ ) (b.4) 

J «^oo ( s (l + s)n) 3 / 2 y ' 

Having in mind the search for j max , 8 d^ ax = max, in deriving (|B .3|> we expanded /(j, 9) as /(j, 6*) w a9 2 + b9 4 . Indeed, 
solving d s d^ /dj = we get 



^«Kl + .),-» _l _ (1 + 2^ (B.5) 

A comparison with the exact results for the degeneracies for s — \ and s = 1 shows an excellent agrement with the 
asymptotics for j < y/n. In both cases numerics show that the asymptotic estimate of j max given above differs from 
its exact value by a factor of \/2/3 and that the maximal degeneracies are of the same order of magnitude. 

APPENDIX C: ENTANGLEMENT CALCULATIONS 

Since Eq. (|A.4fl gives the entanglement of formation of p, it is necessary to find the multiplicities of the spin-j 
subspaces. For s = h the result is 

where instead of the subspace labels A, B we write the number of qubits. The multiplicities satisfy the following 
normalization conditions: 

df k) df n - 2k) = N, £ df k) = , d f k) W + !) = 22 "- ( C - 2 ) 
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We start the entanglement evaluation by considering two equal subsystems of n qubits each. In this case d- : = dj, 



so 

1/2 

X 



E(p\n:n) = ±J2d 2 j log(2j + l). (C.3) 
3=0 

In the leading order 

E(p\n : n) ~ log(2jW x 2 ) + 1) • 1, (C.4) 
where the coefficient c^ 2 ^ reaches its maximal value at 

i=£^V^+2-l, (C.5) 

so 

rc) ~ |logn. (C.6) 

Numerical simulations show that E{p\n : n) — \ logn « 0.0183 

In a generic case of 2k : 2n ~ 2k splitting for 1 <C k < n/2 the sum in Eq. (|A.4|I can be evaluated similarly to 
Eq. (|C.4ll . In this case 



,(„,*) „ 1 / 2 + 3n + 4k(n=k) 
W ~ 2 y n + 1 



« « o \/ — ^ - 1, (C7) 



which reduces to the result of Eq. I|C.5|) for the n : n splitting. It is also interesting to note the fraction of the 
unentangled states in the decomposition of p 

s ( 2fe) _ (C 8) 

goes down with k, 1 -C k < n/2, 

3/2 



,(2fe) 



n 



tt \ k(n — k) 



(C.9) 



The quantum mutual information was introduced in Sec. 2.2. The reduced density matrix of a smaller block for 
2k :2n — 2k splitting is 

k < 

ppriEEfe. ( c - 10 ) 



N ■ 

j=0 a d = l 



and it entropy is 



3=0 

In particular, for the n : n decomposition 



n/2 d 2 AT n/2 d 2 

5(p (n) ) = 25G>) + J] ^ J, = E( P ) +logiV-E# lo S^- ( C - 12 ) 

i=o J ' 3=0 
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Again we approximate 

n/2 j2 n/2 ^ 2 

E # lo s d i ~ E if lo § d *»- = lo s > ( c - 13 ) 

where log dj majt ~ n log 2 — log n, so 

5(p (n) ) - £J(p|n : n) + n log 2 - ± log n + . . . - n log 2. (C. 14) 

As a result, 

I p (n : n) = 2S(p (n) ) - S(p) ~ 3£(p|n : n). (C.15) 

Using the results of Appendix B it is possible to show that the above relation is true for any fundamental spin s. 
Indeed, from Eq. (|B.4fl it follows that = dimH satisfies 

q 

logiV (s) =log s 4 (n) (0)~nlog(l + 2fc)--logn, (C.16) 
J 2 

from (|B.5(1 and 1C.13(I the entanglement between the halves is 

£(p)~ilogn, (C.17) 

and the entropy of one-half of the horizon state is 

77 

S(p n/2 )~-log(l + 2fc)-logri. (C.18) 



As a result, 

3 



/ pS (|:|)~-logn. (C.19) 



APPENDIX D: OVERVIEW OF THE REPRESENTATIONS OF THE PERMUTATION GROUP 

We list some of the basic properties of the representations of S n . An exhaustive discussion can be found, e.g., in 
[l3l FLU 122^ . Irreducible representations of a permutation group are labeled by the Young tableaux, and those that 
correspond to the irreps cr„j consist of just two rows. 



1 


2 






n + j 


n + j + 1 




2n 





(D.l) 



The Young tableau of Eq. IjD.lf) is a normal table, where the numbers 1, 2, . . . , In appear in order from left to 
right and from the upper row to the lower row. The general formulas that define the dimensionality of the irreducible 
representations of the permutation group reduce in the case of a n ,j to dj 2 ™^ whose properties are described in Appendix 
C. The dimension of a representation equals to the number of distinct standard tableaux, where number in each row 
appear increasing (not necessarily in strict order) to the right and in each column appear increasing to the bottom. 
The rule for their ordering will be given below. For each partition [v] we enumerate the standard tableaux as {[v], m), 
where m — 1, . . . , dim <j\ v \ . 

The standard orthonormal basis for the irreducible representation of S n that corresponds to the partition [v] (Young- 
Yamanouchi basis) is labeled by the irreducible representations of the groups in the subgroup chain S n D S^-i D 
. . . D S2 to which it belongs. A simple rule describes this construction: from a given Young tableau {[u],m) of n 
objects one obtains another tableau {[v'\,m!) of n — 1 objects by removing the box that contains the number n, etc. 
For example, a particular basis vector of the [31] irrep of S4 is labeled by the following chain of irreps: 
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2 



(D.2) 
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i.e. this irreducible basis vector belongs the irrep [31] of S4, [21] of S3, and [11] of 52. This labeling uniquely 
specifies the basis vectors. They are ordered by the Yamanouchi symbols (r n r n -i . . . ri), where rj is the row number 
of i. For example, the above vector corresponds to the Yamanouchi symbol (1121). Once all the standard tableaux 
of [v] are labelled by their Yamanouchi symbols, they can be linearly ordered, with the vector with the largest 
Yamanouchi symbol being labelled |[z/], 1), the vector with the second largest symbol being \[v], 2), etc. In practice, 
the irreducible vectors are found as the simultaneous eigenvectors of the complete commuting set of operators. In 
the case of a permutation group, this set consists of the 2-cycle class operators for the all subgroups of the chain 
5„ D S„_i D . . . D $2- 

The matrix elements of all operators that act on the irreducible representation space <T[ U ] of S n can be reconstructed 
from the representations of n — 1 generators of S n , which are the 2-cycles permutations (12), (23),. . . , (n — 1, n). Their 
matrix elements in the | [v] , m) basis are obtained from the following set of rules: 

1. If i — 1 and i are in the same row or column of a standard tableau ([v], m), then 



±\[v],m), 



(D.3) 



with the plus sign taken when both numbers are in the same row. 

2. If i — 1 and i are not in the same row or column, then the matrix clement D^, m is defined with the help of the 
axial distance I. It is calculated as follows: starting from the box that contains i — 1 one proceeds by a rectangular 
route, stepping one box each time, until the box with i is reached. Each step right or upward contributes +1 to the 
axial distance, while a step downwards or left decreases s by 1. The matrix elements 



D^ m = ([u],m'\(i-l,i)\M,m) 



(D.4) 



are given by 



m'm 



1/1, 



m = m 



Vl^~T/l, (M,m') = (*-l,l)(H,m) 
otherwise 



(D.5) 



For a given standard tableau (r„r„_i . . . n) and a permutation (i — the axial distance is 

I = q - Cj_i - (n - r-j-i), 



(D.6) 



where Cj is a column number of the object i. For the two-row Young tables the column numbers are easily obtained 
from the corresponding Yamanouchi symbol as 



2 *-£J-=i»j. n = 1 



(D.7) 



For example 
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(34) 
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(34) 




(D.8) 



An alternative method of constructing irreps of S n is to work on its group algebra CS n . For a given Young tableau 
A = ([i/],m) one defines a subgroup P C S n that preserves the rows of A, and the subgroup Q\ that preserves its 
columns. Then the symmetrizer a\ and the antysymmetrizer b\ are defined as 



g eP 



b\ = Y si S n (5).9- 



(D.9) 



These two quantities define the irreducible (or Young) symmetrizer 



s\ = a\b\, 



(D.10) 
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which generates the irreducible representation A that corresponds to the partition [u] by right multiplication on 
CS n . The symmetrizers that are obtained from different standard tableaux of the same partition lead to the distinct 
equivalent irreps. 
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